Optical Properties and Correlation Effects in Na^Co02 
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We have calculated the optical spectra of Naa;Co02 for a;=0.3, 0.5, and 0.7 within the LDA. We 
compare our results to available experimental data and show that the important features and trends 
are reproduced well, but there is a nearly uniform shift of peak positions and poor agreement in 
intensities. We show, through application of a simple model, that these differences can be attributed 
to overhybridization between Co and O orbitals and spin fluctuations which renormalize the band- 
width. Applying the LDA+U procedure shifts the optical peaks further from their experimental 
locations, indicating that this method of incorporating correlation effects is ill-suited for the case 
Na:i,Co02. 
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The layered cobaltate Naa;Co02 is the subject of con- 
siderable interest not only because of a possibly un- 
conventional superconducting stated, but also an un- 
usual phase diagram that encompasses two very dif- 
ferent metallic states, and at least one insulating 
region- , and unusual magnetotransport and thermoelec- 
tric behavior. '^•■^•^ For x < 0.5, the system behaves as a 
paramagnetic metal with Pauli-like susceptibility and su- 
perconducts when hydrated in the narrow region around 
a;=0.3. At a;=0.5, a charge ordered insulating state, 
which is likely anti- ferromagnetic, occurs. When x is 
increased above 0.5, the metallicity returns, but the sus- 
ceptibility becomes Curie- Weiss likej^ and a spin-density 
wave appearsi around a:=0.75. It has been recently 
reported;^ that, at low temperatures, a second insulating 
state may set in at a;=0.25. 

LDA calculations show^ a band structure that evolves 
smoothly with x, in contrast to the sharply different 
regions described above. The experimental phase dia- 
gram is not well reproduced, and a weakly ferromagnetic, 
metallic ground state is found incorrectly for all Na con- 
centrations. Neglect of strong correlation effects by the 
LDA is the most obvious culprit, but their precise role in 
Naa;Co02 is not understood. Circumstantial evidence in- 
dicates that it is substantial, with most estimates of the 
Hubbard U quite large in comparison with the t2g band- 
width. However, the system does not behave characteris- 
tically as a doped Mott-Hubbard insulator, particularly 
around the superconducting composition where suscep- 
tibility measurements indicate a lack of local magnetic 
moments on Co d ions.^^'^"'^ Also important is the ob- 
servation that LDA calculations in systems with strong 
Hubbard correlations routinely underestimate the ten- 
dency to magnetism (the larger the difference, U — I, be- 
tween the Hubbard repulsion and Hund rule coupling, the 
greater the underestimation)^^, and overestimate the an- 
tiferromagnetic superexchange J. In the LDA, Naj;Co02 
is magnetioi^, with antiferromagnetic instabilities that 
are smaller than ferromagnetic. This requires effects be- 
yond LDA to suppress magnetism and casts doubt on 
frustration as a mechanism for suppression of magnetism. 
Because of this, despite the large U/W ratio, spin fluctua- 
tions may offer a more plausible reason for the deviations 
from the LDA. 



In this context, it is important to verify to what ex- 
tent LDA bands correctly describe the one-electron ex- 
citations and fermiology of Naa;Co02. Quantum oscilla- 
tion probes would be ideal to investigate the Fermi sur- 
face (FS), particularly the predicted small pocketsi^ti^, 
but so far no such studies have been reported. Several 
photoemission reports are availabl o^^i^^ , but because the 
electronic structure is expected to be sensitive to both 
surface states and surface quality, generalization of these 
results to the details of bulk electronic properties is ques- 
tionable. Infrared optics has a penetration depth of ~ 
c/ujp (where LUpi is the plasma frequency), much longer 
than that of photoemission (~ 10 A), and may therefore 
provide the most reliable probe of the electronic structure 
available to date. Several of these experiments have been 
reportedSiiii24i2iS£. However, since optical absorption is 
an integrated property, involving a variety of interband 
transitions, it is usually difficult, if even possible, to in- 
terpret the results in terms of the electronic structure 
without full first principles calculations. 



Here we report such calculations. By comparing our 
calculated LDA spectra with optical experiments re- 
ported for various Na contents, we are able to assign 
the three main peaks to specific interband transitions 
and thereby examine in detail the dependence of these 
peak positions and strengths on bandwidths and ener- 
gies. The peak shifts one expects to see in Mott-Hubbard 
type compounds are not realized in this system, even for 
the insulating compounds. Furthermore, application of 
the LDA-I-U methodology exaggerates the systematic dis- 
crepancies between the LDA and measured spectra. We 
employ a simple model to illustrate that overestimation 
of Co-0 hybridization due to overly extended d-orbitals 
would produce precisely the kind of peak shifts we see in 
our LDA spectra. Spin fluctuations predicted early onSi 
and later observediSS. at some values of x, could also af- 
fect the spectra by renormalizing the bandwidth of the 
t2g complex. This is also consistent with the differences 
we observe in our transitions compared with experiment. 
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I. CALCULATED OPTICAL CONDUCTIVITY 

Using the Wien21i22, full-potential augmented plane 
wave + local orbital code (APW+lo), we calculated the 
optical conductivity for Naa;Co02 at three Na contents: 
a;=0.3,0.5 and 0.7. To vary the Na content, we used the 
virtual crystal procedure of occuping each site with an 
ion of fractional charge. For both c and a lattice con- 
stants, we used experimental datai^ but relaxed the api- 
cal O height separately at each x. The APW sphere 
radii were the same for all calculations: 2.0 for Na, 1.85 
for Co and 1.55 for O. Rk^ax (essentially a measure of 
the plane-wave cut-off) was set to 7.0, giving a basis set 
of 833 APW's and 60 local orbitals. Our density was well 
converged using 480 fc-points in the irreducible Brillouin 
zone (BZ). 

We first calculated the imaginary part of the dielec- 
tric constant, f" {(jj) in the random phase approximation 
(RPA) using the standard code of the Wien2k package. 
The real part of the optical conductivity, a' (w) , can then 
be obtained by (j'{u}) ~ |^e"(u;). The resulting spec- 
tra exhibit three main peaks, which we label 7, (3, and 
a, according to the notation of Wang et al^S,. We can 
identify the specific interband transitions giving rise to 
these peaks using a band by band decomposition of e"(w). 
Each of 7, and a have distinct origins. The 7-peak 
corresponds to transitions between different bands of the 
metallic complex, the /3-peak to t2g-Eg transitions, 
and the a-peak to transitions from fully occupied 0-p 
states to unoccupied Cg states. Note that the 7 and (3 
peaks exist despite the dipole selection rule which forbids 
d-d transitions, because of O hybridization throughout 
both d-band complexes. The t2g bands are split into an 
aig and two e'g bands by the trigonal crystal fieldi^, and 
all 7 transitions are from a lower e^-like state to an upper 
e'g or flig-like state. The Cg-e^ transitions are wholly due 
to small hole FS's near the K points of the BZ and con- 
tribute to the lowest visible peak in the calculated spec- 
trum (the 7 subpeak near 0.4 eV in Fig. 1). However, 
there are allowed e'g-aig transitions in the same energy 
range that contribute roughly equally to the peak inten- 
sity and therefore, the existence of these hole pockets 
cannot be inferred from the existence of the peak. The 
calculated (3 peak has several distinguishable subpeaks 
corresponding to transitions from different parts of the 
BZ. In FigQ] the peaks of the spectrum for Nao.5Co02 
are shown and representative transitions are indicated 
along high symmetry points of the corresponding band 
structure. 



II. COMPARISON WITH EXPERIMENT 

The three (7, /?, a) conductivity peaks are clearly 
recognizable in all available optical data. The re- 
ported positions and strengths of these peaks, however, 
vary somewhat between experiments, even for identical 
dopings^iiSiiS. For consistency of comparison, we use 
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FIG. 1: top panel The LDA optical conductivity for 
Nao.5Co02 with arrows delineating the three main peaks: 
7 (lined arrows), 13 (solid arrows) and a (dotted arrows) A 
Gaussian broadening of 0.06 eV was used, bottom panel In- 
terband transitions corresponding to peaks in the optical con- 
ductivity. 



mainly the report of Hwang et al^, where spectra for 
several compounds of different Na content are available. 
Overall, the LDA spectra are seen to be in good quali- 
tative, but poor quantitative agreement with experiment 
(see Fig. 12), a situation we will attribute to LDA over hy- 
bridization, and/or spin fluctuations in the next section. 
But first, we discuss the doping dependencies and gen- 
eral shape of the spectra. These are well reproduced by 
calculation. There are two clear trends followed by both 
the LDA and measured (Ti(w) as electron count is in- 
creased (Na added). First, the a and (3 peaks shift away 
from each other and become broader. Second, the 7 peak 
sharpens and moves slightly lower in energy. Both effects 
have their origin in an upward shift of Co-derived bands 
as the Co d band filling is increased. The greater en- 
ergy difference between the 0-p states and the shifted eg 
band determines the a peak position, and a reduction 
in hybridization between Co and O orbitals, now fur- 
ther separated in energy, weakens the t2g-e'g crystal field 
splitting, causing the downward /?-peak shift. Since the 7 
-peak is due to transitions within the t2g band complex, 
its sharpening and shift down are both attributable to 
a slightly narrowed t2g bandwidth caused by decreased 
O mixing. Our calculated in-plane plasma frequency at 
LOp — 1.39 eV , agrees well with values extracted 
from experimen^iSiiS: uj^ = 1.48 eV and LOp = 1.17 eV . 
The c-axis plasma frequency is calculated to be 1.48 eV , 
nearly the same as the in-plane value, reflecting the lack 
of anisotropy at higher Na levels^i. As Na is removed, the 
in-plane plasma frequency and the electronic anisotropy 
increase until, at a;=3, the in-plane Up is 2.86 eV and the 
c-axis LOp is 0.337 eV , approximately an order of magni- 
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tilde apart. The increase in Drude weight and concurrent 
enhancement of anisotropy with electron removal are in 
good agreement with experimental observationsSiii. 

Despite these instances of good agreement, some dra- 
matic differences between the spectra are obvious. The 
experimental peak intensities are almost a factor of two 
smaller and approximately 0.5 eV lower in energy than 
our calculated ones when a Gaussian broadening of 0.06 
eV (this value reproduces all three reported peak widths 
well) is used. Experimentally, the strong dip just be- 
low the a peak deepens from x=0.25 to a;=0.5, but then 
becomes shallow again at x=0.7, whereas the LDA gap 
deepens and widens monotonically with increasing Na 
content. It has been previously suggested both /3 and 
a peaks could be identified with transitions across t2g- 
Cg energy gap and that spin splitting of t2g states could 
account for differences between observed peak positions 
and calculated paramagnetic band positions^^. Though 
no observation of static magnetism has been made for any 
X < 0.75, the LDA is known to predict an FM ground 
state for all x, suggesting that magnetic fluctuations may 
be present in the systemii2i. Direct observations of FM 
spin fluctuations^^ and a spin density wave^ at a;=0.75, 
and c-axis AFM orderingS^ at a;=0.85 further support the 
idea that magnetism may play a role in optics. However, 
since regular optics does not allow for spin-flip transition, 
spin splitting may shift optical transitions only if one of 
the two electronic states involved in an transtion is split 
and the other is not. This is not the case for either 7 
or (3 peaks. There is some effect on the position of the 
a peak, because of oxygen admixture, but it is merely 
0.15 eV even at x—0.3. Thus, the observed shifts of the 
absorption peaks compared to LDA, especially the two 
high-energy peaks, a and cannot be due to magnetic 
effects. 



III. CORRELATIONS EFFECTS AND BAND 
RENORMALIZATIONS 



The small t2g and eg bandwidths of Naa;Co02 sug- 
gest that correlation effects in Co-d may be impor- 
tant. Our renormalized atom calculations for Hubbard 
U on Co yield U >3.7 eV, while the t2g bandwidth 
W is 1.5 eV, so that the Hubbard ratio (accounting 
for degeneracy) is > 1.5. Spin-unrestricted 

band structure calculations yield a half-metallic ferro- 
magnetic ground state, contrary to experiment. All this 
suggests that conventional LDA calculations should be 
taken with a grain of salt in this compound, and it 
is quite likely that the actual experimental electronic 
structure is considerably renormalized compared to the 
LDA results. Indeed, angle-resolved photoemission spec- 
troscopy (ARPES) measurementsi^ii^ have detected a 
heavy quasiparticle band near the Fermi energy with a to- 
tal width of ~ 70 meV, much smaller than the LDA band- 
width. This is an interesting finding, but surface sensitiv- 
ity should be kept in mind. Indeed, differences between 
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FIG. 2: (color online) A comparison of experimental (top 
panel) and calculated (lower panel) optical spectra. Note that 
the lowest calculated Na concentration is 0.3, slightly different 
from the experimental value of 0.25. 



surface electronic structure, as measured accurately by 
ARPES, and bulk electronic structure have been estab- 
lished in other oxides, including cuprate^^, ruthenateSi, 
and manganite^SiSS compounds. Optical experiments, 
with their longer penetration depth have the potential 
to tell us more about the specific shortcomings of LDA 
calculations in sodium cobaltate. 

Let us first recall the typical inadequacies of LDA cal- 
culations in systems with localized electrons. The best 
known case is exemplified by Mott-Hubbard insulators. 
Here the missing physics is mainly the Coulomb repul- 
sion between electrons localized at the same site. Typical 
examples are /—electron compounds, high-Tc cuprates, 
or 3d oxides. This group is characterized by (1) un- 
derestimation of propensity toward magnetism (up to 
the level of total loss of a magnetic ground state, as in 
cuprates) (2) underestimation of band gaps between oc- 
cupied and empty bands and (3) overestimation of the su- 
perexchange antiferromagnetic interaction (this increases 
with hopping and decreases with the band gap; the for- 
mer is overestimated and the latter underestimated in 
LDA). As a leading correction to the LDA, the LDA-I-U 
method has been very successful in the prototypical sys- 
tems listed above. For these, the LDA+U method repro- 
duces the missing Mott-Hubbard effects, such as "sink- 
ing" of the occupied bands relative to the unoccupied 
bands (with corresponding shifts of interband transitions 
to higher energies). 

Inspection of the LDA bandstructure of Na2;Co02 in 
conjunction with experimental data indicates that this 
physics is not applicable here. The tendency toward mag- 
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netism is overestimated, and so is the band gap between 
the and Cg bands. The LDA ground state is ferro- 
magnetic, and not antiferromagnetic. Not surprisingly, 
LDA+U calculations22i^ only worsen the situation. The 
t2g — eg distance, already 0.5 eV too large, increases fur- 
ther, and the tendency towards magnetism becomes even 
stronger. Additionally, for more subtle reasons, the two 
e'g subbands of the t2g manifold are shifted down with 
respect to the aig band. This in turn shifts the 7 peak, 
already too high in energy, up even further. We con- 
clude that the LDA + U approach is not appropriate for 
NaxCo02- 

Another possible manifestation of electron-electron 
correlations we is found in '^He and in CrO^^. Here col- 
lective excitations (magnons or paramagnons) play the 
same role as phonons in the sense of "dressing" quasi- 
particles and increasing their massS. As a result, the 
optical spectra more or less proportionally squeeze to- 
ward the low frequency. Softening of magnetic excita- 
tions near a quantum critical point leads to large spin 
fluctuations and suppression (full or partial) of magnetic 
ordering'^^''^^. A comparison of the experimental and cal- 
culated properties of Naa;Co02 indicates the presence 
of such effects. Thus, quantum criticality is likely to 
be an important reason for deviation from the LDA in 
NaxCo02- 

Finally, one should not forget about a very prosaic 
shortcoming of LDA: the fact that it includes a spuri- 
ous self-interaction which leads to overextended 3d or- 
bitals, and thus overestimated hybridization with ligand 
orbitals. This effect becomes stronger for more localized 
d-orbitals. For the narrow d-bands of Naa;Co02, this will 
be operative to some extent. 

We shall now try to assess qualitatively the rami- 
fication of these two effects. The first, "dressing" of 
one-electron excitations is qualitatively similar to the 
electron-phonon coupling induced renormalization, ex- 
cept it occurs in a larger energy range. It is reasonable 
to expect it to affect the whole t2g band or a large part 
of it, leading to overall narrowing of all three t2g bands. 
This, in turn should shift both components of the 7 tran- 
sition, the e'g — e'g one and the e'g — aig one to lower en- 
ergies, without much change in their intensity (the extra 
spectral weight is transferred to high energies, as in the 
electron-phonon coupling). It is unlikely that, as some- 
times assumed, the top of the e'g band will be shifted 
down with respect to the aig band. First, dressing of the 
quasiparticles makes them heavier, but normally does not 
shift different states near the Fermi energy with respect 
to each other. Second, this would not only eliminate the 
e'g — e'g transitions, but would also shift the e'g — aig to 
higher energy, opposite to what is required by the experi- 
ment. It is also unlikely that such "dressing" would affect 
the unoccupied eg band, which is too far away from the 
Fermi level. 

Let us now estimate the possible effects of contraction 
of the Co d -orbitals. To this end we use the simplest pos- 
sible tight binding model, namely one where only nearest 



neighbor Co-0 pda and pdw hoppings (we will these t"^ 
and t^ respectively) are included, and the energy of the 
O p level is taken to be sufficently far removed from the 
Co d level to be integrated out. Both assumptions are 
extremely crude, for direct Co-Co and especially 0-0 
overlaps are not small, nor is the p~ d energy separation 
small, but the model nonetheless provides a useful qual- 
itative framework. The resulting 5x5 Hamiltonian has a 
2x2 Cg subblock with diagonal elements proportional to 
t^. Because there is no path connecting an eg orbital on 
one site to any neighboring eg orbital via oxygen, there is 
no dispersion in this subblock. The 3x3 t2g subblock has 
noticeable dispersion controlled by tj, with the overall 
upward shift of the same order. Finally there is a disper- 
sive off-diagonal subblock, corresponding to 0-assisted 
hopping between neighboring t2g and eg orbitals. The 
scale of this block is set by the product ta-tj^- Note that 
in a cubic structure, such t2g-eg hopping through an in- 
termediary O is impossible. 

We now make a further simplifying assumption, that 
tg. ^ i^. We can then neglect the ligand field on the t2g 
orbitals (proportional to i^,), leaving only the ligand field 
on the eg orbitals (proportional to t^). Since the /? peak 
is due to t2g — eg transitions, its position reflects a crys- 
tal field splitting of the Co d-states which results from a 
combination of electrostatic and ligand field effects. Us- 
ing a linearized muffin-tin orbital (LMTO) calculation 
with the Co-0 hybridization suppressed, we found that 
the crystal field splitting reduces from « 3.4 eV to « 1 
eV. That is, the electrostatic crystal field is 1 eV, and 
the ligand field is 2.4 eV. The latter is probably overesti- 
mated because the LDA overhybridizes the Co d orbitals. 
Sp ecifically , if to- were 12% smaller than it's LDA value 
(v/2.4/1.9 = 1.12), the the /3 peak would shift down by 
0.5 eV, in agreement with the experiment. Furthermore, 
the intensity of this peak would also be reduced, because 
the optical matrix elements in a two-band tight binding 
model scale with the off-diagonal elements of the Hamil- 
tonian (this follows from the tight-binding defininition 
of dipole matrix elements as {Ai\dH/d'k\A2) , where A 
are the eigenvectors), and these scale as to-^ir- This is in 
accord with overestimation of the /? peak intensity and 
lends further support to our assumption that the LDA 
overestimates the pda hopping amplitudes in Naa;Co02, 
probably by some 10%. 



IV. CONCLUSIONS. 

We have calculated the LDA optical conductivity spec- 
tra of Naa;Co02 and compared it with experimental data 
for several different Na contents. There is good quali- 
tative agreement in terms of the number of peaks and 
their behavior with changing x, but exact peak heights 
and positions are not well-reproduced. Though the un- 
derlying reasons for this are likely correlation effects ne- 
glected by the LDA, we argue that using LDA-I-U as a 
remedy is inappropriate for this compound. Spin flue- 
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tuation driven renormalization of the t2g band complex 
and overextension of Co d-orbitals are shown to affect the 
optical spectrum in a manner consistent with the discrep- 
ancies between calculation and experiment. These cause 
narrowing of the t2g band and reduction of crystal field 
splitting respectively, bringing LDA optical peaks into 
good agreement with experimental reports. 
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